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Abstract 

Recent theoretical work has not led to a consensus regarding the nature of the low-energy El strength in the Ca 
isotopes, for which high-resolution (7, 7') data exist. Here we revisit this problem using the first-order quasi-particle 
random-phase approximation (QRPA) and different interactions. First we examine all even Ca isotopes with A'^ = 
14 — 40. All isotopes are predicted to undergo dipole transitions at low energy, of large and comparable isoscalar 
strength but of varying El strength. Provided a moderate and uniform energetic shift is introduced to the results, 
QRPA with the Gogny DIS interaction is able to account for the (7, 7') data, because, up to A' = 28, it yields a rather 
pure isoscalar oscillation. A neutron-skin oscillation is anticipated for A'^ > 30. This contradicts existing predictions 
that Ca develop a neutron-skin mode. Which theoretical result is correct cannot be resolved conclusively using 
the available data. We propose that alpha-scattering, possibly followed by an electroexcitation experiment, could 
resolve the situation and thereby help to improve the different models aspiring to describe reliably the low-energy 
dipole strength of nuclei. 
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1. Introduction 

Recent theoretical work [1-3] has led to a par- 
tial understanding, at best, of the nature of the 
low-energy El strength observed in Ca isotopes [1]. 
Effects beyond first-order random-phase approx- 
imation (RPA), or quasi-particle RPA (QRPA), 
seem necessary to reproduce certain details of the 
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strength distributions, but contradictory reports 
exist as to whether a collective transition is present 
at all in ^**Ca [2,3]. Here we show that QRPA with a 
conventional force, Gogny DIS, also reproduces the 
observed systcmatics and provides a new physical 
interpretation for it. We then propose a strategy to 
resolve which theoretical scenario is correct. 

In the following we avoid the loaded term "pygmy 
dipole" strength, or resonance. We instead refer to 
electric-dipole transitions below the giant dipole 
resonance (GDR) simply as low-energy dipole 
states, or LEDs, regardless of their possible nature, 
and to their strength, especially El, as low-energy 
dipole strength, or LED strength. The acronym IS- 
LED stands for isoscalar low-energy dipole states 
or strength. To refer specifically to an oscillation of 
a neutron skin against a rather inert core, we use 
the term "neutron-skin mode" , or NSM. 
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Let us first recapitulate the current situation. 
High-resolution (7,7') experiments have measured 
the El strength at excitation energies below 10 
MeV in the three isotopes, *°'*^'^^Ca [1]: rather 
little in '**'Ca, and roughly 10 times more in ''^Ca, 
■^^Ca. That there is comparable and even some- 
what less strength in ^^Ca than in ^"'Ca is viewed 
as a puzzle, because it contradicts the simplistic 
expectation, founded on early calculations [4], that 
LED strength should increase with neutron/proton 
asymmetry. The above reasoning disregards the fact 
that the strongest LED state in the N = Z nucleus 
'^''Ca carries too much El strength compared to 
other nuclei of similar mass [5]. Even so, RPA tends 
to predict too much El strength for this state [6]. 
For a useful comparison between theory and exper- 
iment one must therefore consider ''"Ca to be part 
of the systematics, not just a core with negligible 
LED strength. 

We point out that Ca isotopes cannot be consid- 
ered too light to develop collective LEDs. Precisely 
such a mode is observed in ^°Ca (and even in ^^O) 
- see Ref . [6] , where a compelling case was made for 
its coherent nature. If '**Ca does not develop a col- 
lective LED mode at all, its mass is not the underly- 
ing reason. Whether or not it is too light to develop 
a NSM [3] , is a different question. We stress that a 
NSM is not the only collective dipole vibration that 
could develop below the GDR and therefore not the 
only conceivable collective mechanism for generat- 
ing LED strength. This becomes obvious in the (;asc! 
of ''°Ca [6] and was stressed also in Ref. [7], devoted 
to the toroidal dipole mode, and implied in Ref. [8] . 

In general, the first-order RPA and QRPA can- 
not describe precisely the LED strength observed in 
Ca isotopes. In self-consistent (Q)RPA calculations, 
the lowest dipole states tend to lie; too iiigh in exci- 
tation energy, or the systematics tends to be wrong: 
El strength is predicted to rise almost linearly with 
neutron number N [4], apparently due to the de- 
velopment of a neutron skin. Moreover, as already 
mentioned, the LED strength of ^°Ca is often over- 
estimated. The so-called Extended Theory of Finite 
Fermi Systems (ETFFS), which goes beyond RPA, 
is actually the only model which has reproduced the 
properties of LED strength observed in *'^Ca, ^^Ca 
and ^^Ca quite accurately [2]. 

Configurations beyond first-order (Q)RPA can 
affect the LED strength by 1) shifting strength 
to lower energies, where LED strength is exper- 
imentally observed, as first suggested in [9]; 2) 
introducing fragmentation and possible quench- 



ing of strength; 3) introducing additional states of 
two-phonon character, not described by RPA; 4) 
redistributing strength so that much LED strength 
remains above the measurement endpoint of 10 
MeV. All the above mechanisms have been touched 
upon in Refs. [1-3]. A strong statement in Ref. [2] 
was indeed that in the case of ^*Ca phonon cou- 
pling shifts a strong LED state resembling a NSM 
to energies higher than 10 MeV and thus the exper- 
iment missed it. This was proposed as an explana- 
tion of the observed systematics. Notably, an early 
Coulomb-excitation experiment found no evidence 
of a NSM mode below 12 MeV in ^^Ca [10]. 

An important concern remains then that RPA- 
based models predict collective LED vibrations for 
all Ca isotopes, of NSM type or other, which have 
not been identified except for ''"Ca. Obviously, RPA, 
like any model, cannot be expected to account for all 
possible collective nuclear states. The converse is not 
true, however: A collective RPA vibration that re- 
mains incomprehensibly elusive is not a minor prob- 
lem. 

At this point we note that, if mechanism 1) above 
is a dominant one, then a reasonable description of 
LED strength should be possible within (Q)RPA, 
once a uniform energetic shift is introduced. We will 
make such an attempt in this work. 

We may now give an outline of the present work. 
First we present results with different effective in- 
teractions, Gogny and UCOM-based, for the IS and 
El strength distributions of even Ca isotopes with 
N = 14 — 40 and discuss their general features. 
We then make an attempt to describe the exper- 
imental systematics within QRPA, by introducing 
an energetic shift as already mentioned. We find it 
successful when the Gogny DIS interaction is used. 
The reason is that, for TV up to 28, a predomi- 
nantly isoscalar state is predicted, as opposed to a 
NSM. Such a scenario contradicts the ETFFS pre- 
dictions [2]. We suggest that an electroexcitation 
experiment on '**Ca, preceded by alpha scattering, 
could resolve the situation. 

2. Theory 

We employ the self-consistent Hartree-Fock- 
Bogolyubov-Quasiparticle-RPA (HFB-QRPA) , 
which reduces to the self-consistent Hartrec!-Fock- 
RPA for closed-shell nuclei. The ground-state prob- 
lem is solved within a single-particle basis spanning 
15 harmonic-oscillator shells. The same effective in- 
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teraction is used to construct the QRPA equations, 
solved within the HFB basis. For details see [11]. 
We consider a two-body Hamiltonian of the form 



H ^T + V^^ + Vcoui + Vp, 



(1) 



where T is the intrinsic kinetic energy, Vnn a two- 
body nuclear interaction, Vcoui the Coulomb inter- 
action between protons and 



Vp = ts{l + xs)5{f)p"{R) 



(2) 



is a density-dependent contact interaction (f the 
relative and R the center-of-mass position vec- 
tor of the interacting miclcon pair). We employ 
the phenomenological Gogny DIS [12] parameter- 
ization and a unitarily-transformed AV18 realis- 
tic potential, supplemented with a phenomeno- 
logical three-body contact term, which we label 
UCOM(SRG)s,53N- See Refs. [6,13] for informa- 
tion on the latter. For comparison we will also use 
Second-RPA (SRPA) for '^"^'^^Cn, with the full cou- 
pling in the 2p2h space, and with the pure two-body 
UCOM-transformed AV18 potential, UCOMvar 
(UCOM-SRPA). This model is rather well suited 
for the description of the GDR region [14]. Here we 
introduce an energy cut-off of 140 MeV in the 2p2h 
space, so that the IS-LED state of '*°Ca appears at 
approximately 7 MeV, i.e., close to its measured 
value. 

The ISD response is determined by the transition 
matrix elements of the operator 



OiSD = ^e{r 



(3) 



and the electromagnetic response by 



Oei = ^ ^r„Fi„(0„) - 5^rpyi„(np) (4) 



p=l 



in an obvious notation, where the subscripts p and n 
refer to protons and neutrons, respectively. We cal- 
culate the excitation strength, B{El'f). The above 
operators include corrections to explicitly restore 
translational invariance. Within our self-consistent 
(Q)RPA calculations, we obtain practically the same 
values of strength if we use the uncorrected forms 
of these operators, except of course for the spurious 
state, which appears at practically zero energy. In 
SRPA spurious admixtures are unavoidable. 

Electroexcitation cross sections are calculated by 
using the proton transition density, 6pp{r). For the 



longitudinal form factors in plane-wave Born ap- 
proximation (PWBA) we use the convention 



OO 

f 



r^dr. 



(5) 



Within the distorted-wave Born approximation 
(DWBA), the cross section divided by the Mott 
cross section takes the place of the form factor 
squared. 

3. Results 

In fig. 1 we show the ISD and El strength dis- 
tributions of the even Ca isotopes calculated using 
the Gogny DIS and the UCOM(SRG)s.53N interac- 
tions. For ^"^^^C'A, UCOM-SRPA results are shown 
as well. The patterns and trends with respect to neu- 
tron number become clear in this mode of presen- 
tation, where the amount of strength at each value 
of excitation energy is represented by the area of a 
circle. A different scale is used for the IS (discs) and 
El (open circles) strength. 

The theoretical peaks of the GDR are visible as 
the larger open circles in fig. 1. The GDR peaks of 
the A'' = 20 — 28 isotopes are found experimentally 
at an excitation energy of, roughly, 18 — 21 MeV. All 
three models reproduce this fact rather well. In all 
cases we find strong ISD states at the lower end of the 
spectrum. The same observation is made in a previ- 
ous QRPA study using a different SRG-based inter- 
action [11] and is in agreement with existing results 
using Skyrme functionals [15,16]. These results are 
not surprising, because IS-LED states below thresh- 
old exhausting 4-14% of the IS EWSR have been 
observed in various nuclei since years [17]. For ^°Ca, 
the lowest-lying RPA state has been identified re- 
cently as the isospin- forbidden El transition [6], ex- 
perimentally observed at about 7 MeV [18,19]. No- 
tably, the IS strength of the lowest strongly IS state 
of all isotopes and for both interactions is of the same 
order of magnitude, but the same does not hold for 
the El strength of those states. 

Let us examine the three sets of results separately, 
starting with QRPA and the Gogny DIS interac- 
tion. We first point out that the IS-LED strength 
is split in the neutron-deficient isotopes and even 
more fragmented for > 30. For the lowest-lying 
peaks in the N = 18, 20, 22, 24, 26, 28, 30 isotopes we 
find, respectively, 590 (2 peaks below 10 MeV), 832, 
860, 868, 855, 845, 1731 (two lowest peaks) [fm^^] of 
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Fig. 1. Isoscalar dipole (full circles) and El (open cir- 
cles) strength in the Ca isotopes vs. neutron number 
N and excitation energy Ex, for the Gogny DIS and 
the UCOM(SRG)s ,53N interactions, within QRPA, and for 
UCOM-SRPA (up to 22 MeV only). The areas of the circles 
are proportional to the strength. A different scale is used for 
the isoscalar and El strength. The shaded area corresponds 
to the energy region 13.25 - 14.25 MeV (cf. fig.2). 



IS Strength and 26.7 (two peaks), 4.5, 1.9, 1.6, 3.0, 
6.0, 315 (two peaks) [x lO-^e^fm^j El strength. The 
IS-LED strength doubles between "^Ca and ^"^Ca 
and remains large for the heavier isotopes. The El 
strength increases fiftyfold from ^^Ca to ^°Ca and 
remains in that order of magnitude in the heav- 
ier isotopes. At = 30 and beyond the energy 
of the lowest peak decreases. The El strength (but 
not the IS strength) reaches a minimum at A'^ = 
24 and not at A^ = 20. A similar counter-intuitive 
result appears in a recent study of shell effects in 
the isotopic behaviour of LED strength [20]. The 
transition densities of these states reveal that they 
resemble a neutron-skin oscillation only for A^ > 
30, not the lighter N > Z isotopes. The NSM is 
of mixed isospin character and appears strong in 
both the IS and El channels. The stronger IS-LED 
states of the lighter N > Z isotopes show a locally 
isoscalar character instead, resembling that of ""^Ca, 
and therefore carry little El strength. In the case 
of '*^Ca, for example, eight out of the 18 possible 
Ihu) ph configurations (four of them proton config- 
urations) contribute at least 1/18 each to the norm 
of the IS-LED state, making it quite a collective 
state [21]. Many ph configurations contribute coher- 
ently to its IS strength. Dipole states in ''^Ca with 
more-pronounced neutron transition densities are 
found at somewhat higher energies, but none can be 
viewed as an oscillation of a neutron skin against an 
inert core. Remarkably, the hole states in ^^Ca are 
found very similar for protons and neutrons, a situa- 
tion that is not realized in '^"Ca. States dominated by 
proton transitions are found in the neutron-deficient 
isotopes. A variety of states with locally isoscalar 
transition densities, where both nucleon species con- 
tribute at the surface with non-negligible tails, are 
present regardless of A^. 

The results for the UCOM(SRG)s,53N interaction 
are more in line with earlier predictions [4]: a min- 
imum in El strength is reached for "^^Ca, and a 
monotonic behaviour is observed on either side of 
A^ = 20. At both A = 22 and A^ = 30 the posi- 
tion of the lowest peak shows a kink, pointing to 
structural changes of this state. Here a neutron-skin 
mode, characterized by a translational-type 6pp{r) 
and an extended Spn{r) with a node, is predicted al- 
ready for the lightest N > Z isotopes. In the A = 
20, 24, 28 isotopes, however, these states carry too 
much El strength: they are about one order of mag- 
nitude stronger than the experimentaly observed 
transitions below 10 MeV [1] . The same holds for the 
results of Ref . [4] . Therefore the NSM scenario fails 
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to explain the observed systematics, at the (Q)RPA 
leveL 

Finahy we look at the UCOM-SRPA resuhs. As 
with QRPA and the UCOM(SRG)s,53n interaction, 
we find that the lowest strongly IS state alone carries 
too much El strength to be compatible with the 
data, namely 0.021 and 0.4 e^fm^ for "OQa and ^^Ca, 
respectively. In the case of ^®Ca, this state is of NSM 
type and dominated by (yhl'i)^^ transitions. Again 
the NSM scenario fails to describe the data. 

Within RPA with the UCOM(SRG)s,53N inter- 
action as well as within the UCOM-SRPA model, 
we find that other IS states of ^®Ca, at somewhat 
higher energies but still below the GDR, are of lo- 
cally isoscalar character with a node in the transi- 
tion densities of both nucleon species. 

From the present results we conclude that all iso- 
topes develop strong IS-LED states, but the nature 
of these states changes as iV increases from resem- 
bling the oscillation of a proton skin, to a rather pure 
isoscalar oscillation, to a neutron-skin oscillation. 
The neutron number at which the latter transition 
takes place depends on the interaction. The energy 
and fragmentation patterns also change. An impor- 
tant observation is that El strength is found in the 
whole energy range from the lowest-lying state up to 
the GDR. Especially for the heavier isotopes, a clus- 
tering of strength in two or more energetic regions 
is apparent. Judging from the variety of the calcu- 
lated transition densities, we expect that the mix- 
ing and fragmentation of different possible vibra- 
tional modes gives rise to this result: oscillations of a 
strongly decoupled neutron (or proton) skin against 
a core, or oscillations of a mixed layer of protons and 
neutrons against a core, these being of toroidal type 
or of partially compressional type - see, for example, 
[8,6,7,22-24] for studies of various possible dipole 
modes with different approaches. In general, the po- 
sition of the torus or compression point may vary 
with energy and TV. GDR admixtures and Ipl/i tran- 
sitions are also possible. A detailed study of all these 
states goes beyond the scope of this work. Neverthe- 
less, we take the opportunity to point out again that 
El strength alone is not a good indicator of collec- 
tivity, in particular the lack thereof. A collective IS 
dipole state could in principle carry less El strength 
than a Ipl/i transition [6]. It may be telling that 
LED strength is often overestimated, even when it is 
interpreted as due to non-collective transitions [25]. 

We now return to our main topic, which is the 
LED strength of 40^44,48^^^^ ^gj^ whether the 
scenario predicted by the use of the Gogny DIS 
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Fig. 2. Summed electric dipole strength for the three isotopes 
up to the indicated excitation energies: experimental data [1], 
ETFFS results [2] and present QRPA results with the Gogny 
DIS interaction. For the Gogny DIS model the El strength 
of only the lowest IS-LED state is also given. 



interaction is compatible with the (7, 7') data. A 
revealing comparison is made in fig. 2. Here we show 
the summed El strength observed below 10 MeV, 
the existing prediction within ETFFS, and our new 
sets of results: The El strength of the lowest dipole 
state in each isotope within (Q)RPA using the 
Gogny DIS interaction; the summed strength below 
13.25 MeV, 13.75 MeV, and 14.25 MeV within the 
same model. The lowest states are too weak with 
respect to experiment. A similar observation was 
made in ref. [16], where the Skyrme SkM* interac- 
tion was used, and in ref. [11]. The energy cut-ofF 
strongly influences the quality of the comparison. 
In fact, minor uncertainties in the data and the the- 
oretical results can likewise affect the systematics. 
Therefore an attempt to reproduce the data very 
precisely in this region appears rather meaningless. 
Nevertheless, the measurements are rather well re- 
produced by QRPA for the higher cutoffs, indicated 
in fig. 1 by the shaded area. This generic observa- 
tion remains true if we change certain ingredients 
of our calculation such as the harmonic-oscillator 
parameter. That it takes an energy shift by about 
3 — 4 MeV for the QRPA to describe the data, is in 
excellent agreement with the early calculations of 
Ref. [9]. Higher-order configurations are expected 
to shift the strength to lower energies and could 
introduce the fragmentation necessary to describe 
the El values more precisely [2,3,9,26,27]. We con- 
clude that the scenario whereby a strong isoscalar 
mode with little El strength is present in ''^Ca, 
and not a NSM, is very well compatible with the 
measurements. 
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4. Resolving the situation in "^^Ca 

We have seen that a strongly isoscalar state de- 
veloping in *^Ca is as compatible with the exist- 
ing data as a mixcd-isospin NSM appearing above 
10 MeV. We now ask how the correct scenario can 
be resolved. Measurements of the full dipole spec- 
trum of ''^Ca could be useful, because a NSM would 
be visible as a strong El transition close to particle 
threshold. Such information can be obtained from 
relativistic Coulomb excitation in proton scattering 
under extreme forward angles as recently demon- 
strated for the case of ^^^Pb [28] . Data for ^^Ca have 
been taken [29] and arc presently analyzed. But if 
no such strong state is observed, the problem re- 
mains unresolved. We now note that both types of 
mode, namely pure IS or NSM, would carry much 
ISD strenght and therefore would feature promi- 
nently in a-scattering. In particular, {a, a'j) is very 
well suited to detect IS-LED strength below thresh- 
old [17,30-32], but has not been used on Ca isotopes 
other than ^"^Ca [19]. Such an experiment on ^^Ca 
would be most useful in resolving the present sit- 
uation. Inelastic a— scattering above threshold has 
been performed on ^*Ca [33]. A resonance beyond 
16 MeV was found, but no strong state was clearly 
observed in the energy region of the NSM mode as 
predicted by ETFFS. The onset of another concen- 
tration of IS strength below 10 MeV could be in- 
ferred from the data, but not conclusively. 

If a strong isoscalar state is found below 10 MeV, 
it would likely not be of NSA1 type. An elcctroexcita- 
tion experiment could then reveal its nature. Inelas- 
tic electron scattering has been performed on ^^Ca 
by different groups, but the analyses have mostly 
focused on 1+ states and other multipolarities. In 
fig. 3 we show the transition densities of the IS-LED 
state as predicted within RPA by the Gogny DIS in- 
teraction and of the NSM mode as predicted, in this 
case, by UCOM(SRG)s,53n. We also show the corre- 
sponding longitudinal form factors within DWBA. 
The NSM has a large form factor already at low mo- 
mentum transfer, while a completely different result 
is obtained for the other IS-LED state. A minimum 
at low momenta would make this state analogous to 
the IS-LED state of ^°Ca [6]. Significant differences 
between the transverse form factors of the two dif- 
ferent modes are also expected. A more dedicated 
study shall be the subject of upcoming work. 
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Fig. 3. Transition densities and longitudinal electroexcitation 
cross section of the IS-LED state in *®Ca based on the Gogny 
DIS and the UCOM(SRG)s,53n interactions. 

5. Conclusion 

Our QRPA results suggest that all even Ca iso- 
topes with N = 14-40 develop strong IS-LED 
states. IS and El strength is present in the whole en- 
ergy region up to the GDR. The nature of the states 
varies with the neutron number and the energy. Re- 
sults with the Gogny DIS interaction suggest that 
the lowest dipole transition becomes a neutron-skin 
mode only for A'' > 30. This result is compatible 
with the systematics of El strength observed in 
40,44,48(^0^^ but contradicts an earlier interpretation 
of the data [2]. We propose that a— scattering fol- 
lowed by an electroexcitation experiment on '''^Ca 
could resolve the situation conclusively. We take the 
opportunity to stress that, in order to understand 
the nature of the whole low-energy dipole spectrum, 
one must examine not just the El strength, but at 
the same time the response to isoscalar operators, 
in particular toroidal and compressional ones [24]. 
These remarks are relevant for all isotopic chains 
and for problems including the nature of pygmy 
dipole strength in very neutron-rich nuclei, the so- 
called isospin splitting of pygmy dipole strength and 
the origin of El and isoscalar strength throughout 
the region below the GDR. 
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